Optical and spectral properties of quantum domain-walls in the 

generalized Wigner lattice 
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We study the spectral properties of a system of electrons interacting through long-range Coulomb 
potential on a one- dimensional chain. When the interactions dominate over the electronic band- 
width, the charges arrange in an ordered configuration that minimizes the electrostatic energy, 
forming Hubbard's generalized Wigner lattice. In such strong coupling limit, the low energy ex- 
citations are quantum domain-walls that behave as fractionalized charges, and can be bound in 
excitonic pairs. Neglecting higher order excitations, the system properties are well described by an 
effective Hamiltonian in the subspace with one pair of domain-walls, which can be solved exactly. 
The optical conducitivity o(lo) and the spectral function A(k, u) can be calculated analytically, 
and reveal unique features of the unscreened Coulomb interactions that can be directly observed in 
experiments. 



I. INTRODUCTION 

Quantum many-body systems often exhibit complex 
behavior, arising from the strong interactions between 
the individual degrees of freedom. In some cases, the ex- 
otic physical properties can be simply explained in terms 
of suitable renormalized collective excitations. Well- 
known examples are the charge-spin separation in the 
Luttinger liquid with gapless bosonic excitations^ soli- 
tons on a CDW condensate* 2 - and fractional charges in 
quantum Hall systems £ 

Another system that exhibits emerging complex be- 
havior is the generalized Wigner lattice (GWL), de- 
fined as the classical charge pattern that minimizes the 
Coulomb repulsion on a discrete lattice. The GWL was 
introduced by Hubbard 4 to explain the charge ordering 
observed in the TTF-TCNQ organic salts, and has since 
been invoked in several classes of narrow-band quasi- 
one-dimensional compounds£&£ Increasing the electron 
bandwidth smears the classical charge distribution, and 
eventually drives the system towards a small amplitude 
charge density wave 4 

However, in the strongly interacting limit, where the 
bandwidth is much smaller than the electrostatic repul- 
sive energy the quantum charge distribution remains 
very close to the classical configuration. In this limit, the 
low- lying excitations are pairs of domain- walls (kinks and 
anti- kinks), that carry fractional charge 4i2jIS It follows 
that the low energy properties of the system can be deter- 
mined by solving the problem of two interacting domain- 
walls, which is equivalent to the Coulomb problem on a 
one-dimensional chain. Since kinks and anti-kinks have 
opposite charges, the lowest lying excitations of the GWL 
will be bound pairs, followed by a continuum of unbound 
domain-walls. Restricting to the subspace with only one 
kink/anti-kink pair, the quantum melting of the GWL 
was estimated in RefiS as the point where the gap in the 
excitation spectrum vanishes. The result is in agreement 



with previous variational estimates^ as well as with more 
recent exact diagonalization data^ indicating that the 
single pair approximation captures the essential physics 
of the GWL phase. On the other hand, it clearly breaks 
down close to the quantum melting transition, where a 
proliferation of domain-walls is expected. 

In this paper, we take advantage of the exact solu- 
tion of the discrete Coulomb problem given in Refsj n > 12 
to evaluate the optical conductivity of the generalized 
Wigner lattice at simple commensurate fillings n = 1/s. 
The formation of excitons manifests through the emer- 
gence of a series of sharp peaks, followed by a strongly 
asymmetric absorption band due to continuum of scat- 
tering states. Remarkably, both the discrete peaks and 
the asymmetry of the absorption continuum are direct 
consequences of the long-range interactions among the 
electrons in the original model, that disappear when the 
Coulomb interactions are replaced by short-range poten- 
tials. 

An analogous formalism is used to calculate the spec- 
tral function A(k,u>). However, since the addition or re- 
moval of an electron to a GWL at filling n = 1/s is equiv- 
alent to the creation of s domain-walls* 9 the single pair 
approximation only applies to the special case n = 1/2. 
The interactions being repulsive, because the domain- 
walls have equal charge, the low-lying excitations in the 
spectral function are unbound scattering states, while 
anti-bound states appear as a set of dispersive quasi- 
particle peaks above the continuum. 13 

The paper is organized as follows. In Sec. II we in- 
troduce the one-dimensional model for spinless fermions 
with long-range interactions, and solve it in the narrow- 
band regime, restricting to the subspace with one pair 
of domain-walls. In Sections III and IV we use this ana- 
lytical solution to calculate respectively the optical con- 
ductivity and the single particle spectral function. The 
results are briefly discussed in Sec. V in connection with 
existing experimental work. 
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II. MODEL AND SINGLE PAIR 
APPROXIMATION 

We consider the following Hamiltonian for fermions on 
a linear chain, interacting through long-range Coulomb 
forces: 



H 



■*E (4^+1+4^-1 



V ^ ( n i — n ) ( n j — n ) 



■3\ 



Here t is the nearest-neighbor hopping amplitude and 
V sets the energy scale of the long-range repulsion, 
V m =V/\m\. {c\,Ci} are, respectively, the fermionic cre- 
ation and annihilation operators, i%i — c\ci is the occupa- 
tion number at site i and n is the average charge per site. 
The spin degrees of freedom are explicitly neglected, cor- 
responding to the limit of large on-site repulsion U. The 
lattice parameter is set to unity. 

Since we are interested in the strong coupling regime, 
it is useful to start by describing the classical solution, 
obtained for t = 0. At simple commensurate fillings 
n = 1/s, the ground state takes one of the s equiv- 
alent configurations that minimize the electrostatic en- 
ergy, with one electron on every s-th site^ The excited 
configurations of lowest energy correspond to the inclu- 
sion of a domain with a different (shifted) ground state, 
delimited by a kink and an anti-kinki^ The energy cost 
to form a classical domain of finite lenght sd (d is the 
number of perturbed unit cells) can be calculated as 



d oo 



A M = E E (Wi 

p—1 m—p 



Vsm-1 ~ 214, 



(2) 



It defines an effective interaction between the domain- 
walls, and can be evaluated exactly in terms of Digamma 
functions. However, a very good approximation is ob- 
tained by expanding in powers of 1/ sd and keeping only 
the leading term 



A s (d) 



s 2 sd 



(3) 



The above expression can be directly interpreted as the 
Coulomb attraction between two defects of charge ±l/.s 
at distance sd, with £ s = — [l — - cot (-)] the energy 
cost for creating two isolated defects at infinite distance. 
The approximate result Eq.|(5J) is quite accurate even at 
short distances^ Unless otherwise specified, we shall ne- 
glect the higher order terms (dipolar, multipolar interac- 
tions) in the following. 

According to the above discussion, the gap in the ex- 
citation spectrum is given classically by A s (l). Actually, 
this can be viewed as the effective coupling parameter 
of the model, that determines the stability of the classi- 
cal charge pattern against quantum fluctuations, so that 
the melting of the GWL is expected when t > A s (l)£ 
Since the gap scales as ~ V/s 3 ~ Vn 3 , we see that for 
a given value of the hopping integral t, the formation 



of a GWL becomes unfavorable at very low fillings. In 
this case an ordinary Wigner crystal phase is more likely 
realized, where the electron wavefunctions are spread on 
lengthscales much larger than the lattice spacing 4 A sim- 
ilar argument shows that the GWL phase is unfavored 
when moving away from the simple commensurate fill- 
ings n = 1/s. Indeed, at generic rational fillings n = r/s, 
it is easily shown that the gap still scales as ~ V/s 3 (i.e. 
it is governed by the periodicity s of the classical pat- 
tern), which can be arbitrarily small regardless of the 
value of nJ^ In the following, we shall therefore restrict 
to the special fillings n = 1/s, where the physics of the 
GWL is more relevant. 



Let us denote \m,d) the classical state corresponding 
to a kink/antikink pair of length sd, with center of mass 
at m. The total two-body wavefunction in the quantum 
case can be written as 



K.X 



Ljs L/2s 

E E e iK "' 

m — — L/s d— — L/2s 



<p K ,\(d)\m,d) (4) 



where the center of mass of the pair moves freely with 
momentum K, while the radial part obeys the following 
discrete Schrodinger equation-!*^ 



■2t(K)[tp(d+l) + cp(d-l)]- 



V/s 3 



E 



(p{d) = 



(5) 

with the boundary condition <^(0) = 0. The factor 
t(K) = tcos(Ks/2) reflects the fact that at each elec- 
tron hopping event, the center of mass of the defect pair 
moves by is/2 lattice sites. The discrete hydrogen prob- 
lem on the semi-infinite chain admits an exact analytical 
solutio n 11 ' 12 in terms of Gauss-hypergeometric functions 
(or, equivalently, of incomplete beta functions). The ex- 
citation spectrum consists, as in the ordinary hydrogen 
atom, of an infinite series of bound states (kink/anti-kink 
excitons), followed by a continuum of scattering states. 
Observing that the two subspaces of different polarity 
d > and d < are disconnected when only nearest- 
neighbor electron hopping is allowed^ the solution of 
the Coulomb problem Eq.([5]) on the infinite chain is ob- 
tained straightforwardly by constructing combinations of 
even or odd parity of the wavefunctions on the half-chain. 

a. Domain-wall excitons. The bound states have 
energies: 



E K , n = H a -Jl6t(K)* 



V 



s 3 (n + 1) 



(6) 



with n = 0, 1, 2, . . . , oo. Introducing the parameter 
sinh(fc„) = V/[4s 3 t(K)(n + 1)], the normalized radial 
wavefunctions read: 

VK,n(d) = A n d e- k " d 2*1(1 - d, -n\ 2; 1 - e 2fe ") (7) 
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where 2-F1 is one of the Gauss hypergeometric functions 
and A n is a normalization constant: 



A n = 2 sinh(fc n )-\/tanh(fc n ) e 



—nk~ 



(8) 



The exciton radius, defined as r n — d \fK,n{d)\ , 

can also be evaluated^ 



(n + 1) 



2 + cosh(2fc„ 
sinh(2fc„) 



(9) 



In the strong coupling regime, r n w (n + 1) and the 
first bound state (n = 0) is strongly localized on one 
site. At large quantum numbers the radius increases as 
r„ ~ (6s 3 t/V)(n + l) 2 approaching the free scattering 
states. 

b. Scattering states. The energy of the scattering 
states is given by the sum of the energies of the indi- 
vidual domain- walls, with k the relative momentum of 
the two-body system: 

E Ktk = S s - 2t [cos(Ks/2 + k) + cos(Ks/2 - k)} . (10) 

The corresponding wavcfunction is>^ 



= B- 1 (*)de iM 2i'i(l - d, 1 - Mj(fc); 2; 1 



(11) 



where 77(fe) = V r /[2s 3 t(i^) sin/c]. The normalization fac- 
tor is determined by the asymptotic behavior at large d, 
and is given by 

B-\k) = g , ^ (fc ) sin 2 (fc) e-^"W (12) 
L smn[7T?7(fe)J 

in the case of an open chain of length L/2s. 

The excitation spectrum is illustrated in Fig. [1] for 
t/V = 0.05 at fillings n = 1/2 and n = 1/3. The gap in 



n=l/2 
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FIG. 1: Excitation spectrum vs. center of mass momentum K 
in the single pair subspace for t/V = 0.02 at fillings n= 1/2 
(left) and n = 1/3 (right). The full lines correspond to the 
3 lowest bound states, while the dashed lines indicate the 
boundaries of the kink/anti-kink continuum. 



the excitation spectrum is obtained from Eq. ([6]) setting 
K = and n = 0: 



A° pt = S s - to 




(13) 



This value is lower than the classical value A s (l) because 
of the energy gain due to charge derealization. The con- 
dition A° pt — > can be used to locate the melting of the 
GWL as the region where charge defects are expected 
to proliferated For the quarter-filled case s = 2, the esti- 
mate t/V ~ 0.12 is in good agreement with the numerical 
results of Ref 



III. OPTICAL CONDUCTIVITY 

The finite frequency absorption at T — can be ex- 
pressed in terms of the current-current correlation func- 
tion through the Kubo formula: 



1 Im y- \(*kMMgs)\ 2 
Lui m e-' u— E K .x + i0 + 



(14) 



K,\ 



where the sum runs over all eigenstates, and the ground 
state energy is set to 0. For nearest- neighbor hopping, 
the current operator on the lattice is defined as j — 
it^2i( c t c i+i ~ c t+i c i)- The above dimensionless expres- 
sion should be multiplied by a prefactor <tq — e 2 a 2 /Kv to 
restore the proper units, where a is the intersite distance 
on the chain and v is the volume of the unit cell. 

Approximating \tpGs) m the GWL phase at fillings n = 
1/s with the classical state \ipoo) = IL c L|0)i we see that 
the effect of the operator j is to create kink-antikink pairs 
of length d — ±1 at any of the m possible sites. The 
matrix elements with the single pair states of Eq. ((4]) are 
thus given, in terms of the eigenfunctions of the semi- 
infinite chain, as 

KiM5'|**,a>| 2 = t 2 -$K,o I^,a(i)I 2 - (15) 

s 



A. Coulomb potential 

The optical absorption can be divided in two parts a = 
fexc + fsc, which represent respectively the sharp transi- 
tions from the ground state to the excitonic states and an 
absorption band due to transitions to the kink/anti-kink 
continuum. The corresponding energies are obtained 
from Eq.© and Eq.^ by setting K = 0. The exci- 
tonic part is an infinite series of delta function peaks 



Vo,n(l)| 2 <5 (w 



0,n) 



(16) 



where the spectral weights are given by Eqs. ([7]) and 
above, making use of the property 2F1 = 1 at d = 



4 



o.oi 



Setting k(uj) = arccos [—(to — T, s )/4t] and converting the 
sum over states in Eq. (fT4")) into an integral, we obtain 
the following analytical expression for the continuous ab- 
sorption band: 



4s A w 



V k(w) 
2s 3 t sin k(cu) 



V 



(17) 



1 



2sH sin k{uj) 



Note that since we are assuming t < V/s 3 , the denom- 
inator can be set equal to 1 up to exponentially small 
corrections (these become important close to the melting 
of the GWL, where anyway the single pair approximation 
breaks down). 

The optical absorption determined above is illustrated 
in Fig. [21a and agrees remarkably well with the exact 
diagonalization results of Ref^ The main discrepancy is 
a shift in the position of the first excitonic peak, that can 
be entirely ascribed to the the use of the pure Coulomb 
potential in Eq. A transfer of spectral weight 

from the excitons to the kink/anti-kink continuum takes 
place as t/V increases, as shown in Fig. [2Jb, where we 
have plotted the quantities J du>a exc (uj) / J dioa(uj) and 
/ duia sc (uj)/ J dujcr(ui). Note that formulas ([To]) and (fT7| 
obey the following property, that can be checked by di- 
rect evaluation: 



(oj) = / duia(uj)u! — irt 2 / s. 



(18) 



This relation shows that the "average" absorption fre- 
quency is independent on the coupling constant and can 
be used to measure directly the value of the hopping in- 
tegral in an optical absorption experiment. 

In the above derivation it was assumed that the poten- 
tial follows a pure Coulomb law, Eq.Q, which allowed 
us to obtain explicit analytical formulas for the optical 
absorption. However, the discrete Schrodinger equation 
([5]) with the exact potential A s (d) of Eq. ([2]) can be solved 
to arbitrary accuracy with small numerical effort. The 
result is presented in the inset of Fig. Oa. We see that 
the analytical expression ([IT)) agrees quite well with the 
numerical result. Discrepancies arise at very small values 
of t/V , where the excitons become localized on very short 
lengthscales, which is where the potential A s (d) sensibly 
deviates from the Coulomb law. 

Note that an expression similar to Eq. (fT7]) was de- 
rived by Gallinar— in the half-filled case n = 1 for elec- 
trons with spin, in the limit of a large on-site repulsion 
U. In that case the lowest-lying excitations correspond 
to pairs of doubly occupied/empty sitesJ^ The energy 
cost for creating such pairs is set by U, while the long- 
range tail of the Coulomb potential acts to bind the op- 
posite charges together. 25 Similar effects have also been 
discussed in the context of conjugated polymers, in which 
case the energy scale of an electron-hole exciton is set by 
the Peierls dimerization gapi& (see also Refill). 
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FIG. 2: a) Optical conductivity of the generalized Wigner 
lattice in the single pair approximation as given by Eqs. (|16[) 
and (fT7|) at n = 1/2 and t/V = 0.02. A finite broadening 
has been introduced for graphical purposes. The open circles 
are the exact diagonalization data of Ref.— . In the inset, the 
same curve is compared with the result obtained with the full 
potential of Eq.© (dashed lines). 

b) Fraction of the spectral weight carried respectively by the 
excitons (full lines) and the kink/anti-kink continuum (dashed 
line) . 



B. Screened potentials 

The strong asymmetry of the optical lineshape Eq. 
(fT"7f ultimately follows from the long-range nature of the 
Coulomb potential. Indeed, although the energies of the 
scattering states are given by the non-interacting expres- 
sion (|10p , their wavefunctions strongly differ from the free 
plane waves especially at short distances (the scattering 
states are required to be orthogonal to the bound states). 
This is reflected in the optical absorption, which involves 
the wavefunction precisely at d = 1, where the condition 
of orthogonality with the excitonic states is most strin- 
gent. To illustrate this issue, we have considered the case 
of a screened potential V/d — > V exp (—d/l) /d. This po- 
tential is convex and ensures the stability of the classical 
configuration considered here,— which is preferable to the 
case of arbitrarily truncated Coulomb interactions, for 
which the ground state is not always univocally deter- 
mined. 

The results obtained for generic values of the screen- 
ing length are reported in Fig. [3] at filling n = 1/2, for 
t/V = 0.02. Upon reducing the screening length, the 
optical absorption shifts to lower frequencies, the asym- 
metry of the lineshape becomes less marked and the ex- 
citonic peaks progressively disappear (below a critical 
value, no bound states are possible and a exc = 0). In 
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FIG. 3: Optical absorption for a screened potential at n = 1/2 
and t/V = 0.02. The labels indicate the value of the screen- 
ing length £. The dashed line corresponds to the unscreened 
Coulomb potential. 



the limit i — » 0, a truncated nearest-neighbor poten- 
tial is recovered. In that case the effective interaction 
A(d) — Ve^ 1 ^ = V\ is constant for all d and the scat- 
tering states reduce to free-particle states. The corre- 
sponding current-current correlation function has a semi- 
circular shape and the resulting absorption band is 

«-<«)- ElR^P) 1 < 19 > 

To conclude this section, let us stress that the exper- 
imental observation of an asymmetric lineshape of the 
form shown in Figf2]a can be taken as a signature of the 
importance of long-range electron-electron interactions in 
a given material, and calls for theoretical models that go 
beyond the intersite repulsion term usually considered^ 



IV. SPECTRAL FUNCTION AT HALF FILLING 

In a photoemission experiment, an electron is added or 
removed from the system in a sudden process that is de- 
scribed by the single particle spectral function A(k, uS) — 
— -ImG (fc, cj). In the spirit of the strong coupling ap- 
proach presented in the preceding section, adding an elec- 
tron to the classical configuration at n = 1/s creates s 
kinks of charge 1/s. Analogously, removing an electron 
creates s anti- kinks of negative charge —1/s. In the spe- 
cial case s = 2, two defects are created and the problem 
becomes equivalent to the one treated previously, albeit 
with repulsive interactions (the analogy is no longer true 
for larger values of s, where the creation or removal of 
an electron take the system away from the single-pair 
subspace). 

Classically, since the interaction potential is repulsive, 
the state of lowest energy in the system with one particle 
added (removed) corresponds to configurations where the 



two kinks (antikinks) are far apart. The creation energy 
of each kink is S 2 /2 = V/4 (cf. Eq. ©). 

In the quantum case, we can write the 1D-Schr6dinger 
equation for the two-body problem, which is formally 
identical to the one presented in Sec. II, except for a 
change of sign in the interaction potential, which results 
in a symmetrically reflected excitation spectrum: the low 
energy part now corresponds to the domain-wall contin- 
uum, while the discrete levels are moved to higher ener- 
gies. Such anti-bound states form, as for the excitonic 
states encountered previously, due to the long-range na- 
ture of the repulsive potential. 
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FIG. 4: Spectral function A(k,u) for s = 2 and t/V = 0.05. 
The different curves correspond to different values of fc, as 
indicated on the left axis. 

The spectral function for w > 0, which describes the 
inverse photoemission processes, can be calculated as 

A{q,u) = J2\(*k,x\c+\iPgs)\ 2 5(u ~ E K ,x) (20) 

where c+ = J2 m e 1 '"^ and \\&k,\) are the solutions 
of the repulsive Schrodinger equation ([5]) . An analogous 
formula holds for the spectral function at w < 0, with 
A(q, —uj) = A(q,uj). Replacing as usual the ground state 
with the classical configuration at t = enforces the con- 
straint that electrons can only be created on the unoc- 
cupied sites. The resulting expression for the spectral 
function is 

A 

where the sum runs over both the continuous spectrum 
(A = k) and the discrete anti-bound states (A = n), 
whose energies are given by Eqs. (fT0|) and ([6]), that we 
reproduce here (in the latter, the quantity E — Yi s changes 
sign due to the repulsive nature of the interactions) : 

E q , n = \ + ^W^^Q" (22) 
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The results are reported in Fig. |4j The width of the 
continuous spectrum is modulated by q and vanishes at 
q = ±7r. The anti-bound states are clearly observable as 
well defined quasi-particle peaks in A(k, uS), that follow 
the roughly sinusoidal dispersion relations Eq. (j22"|) . 

V. DISCUSSION AND CONCLUSIONS 

In this work, we have examined the spectral properties 
of spinless electrons interacting through the long-range 
Coulomb potential on a one-dimensional lattice, which 
constitutes a minimal model to address the effects of elec- 
tronic correlations in narrow-band solids, away from the 
most studied half-filled case. In the strongly interacting 
regime, the charges form a generalized Wigner lattice, 
whose elementary excitations are fractionally charged 
domain- walls, themselves interacting through long-range 
forces. Taking advantage of the exact solution of the 
Coulomb problem on a one-dimensional chain, we have 
derived an analytical expression for the optical conduc- 
tivity at simple commensurate fillings n = 1/s, as well as 
for the single particle spectral function in the special case 
n = 1/2. Both are in good agreement with the available 
exact diagonalization data at n = l/2i 10 i 13 

The sharp peaks emerging in the optical conductiv- 
ity, signaling the formation of domain-wall excitons, the 
asymmetric lineshape of the absorption continuum, as 



well as the sharp quasi-particle peaks in the single par- 
ticle spectral function, are all robust features that follow 
from the long-range nature of the Coulomb potential, and 
that are lost in models with strongly screened or trun- 
cated interactions. These features should be clearly ob- 
servable in experiments, allowing to address the relevance 
of long-range interactions in the charge-ordered insulat- 
ing phases of narrow-band one-dimensional systems. 

As an example, the mid-infrared optical absorption 
band recently measured^ in the organic quarter-filled 
(ri = 1/2) salt (DI-DCNQI)2Ag — a prototypical system 
exhibiting Wigner crystal ordering^ — resembles very 
closely the theoretical spectra obtained from Eq. (TIB"]) 
and (jTTJ) , and a satisfactory fit can be obtained with phys- 
ically sound parameters (V = l.leV and t = 80mey). 
Absorption shapes similar to the ones presented in this 
work have also been reported in other classes of organic 
compounds, such as DBTSF-TCNQF 4 , TTF-TCNQ 21 
and Cs2(TCNQ) 3) 22 corresponding respectively ton — 1, 
n ~ 1/2 and n = 1/3, whose detailed analysis is post- 
poned to future work. 
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